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Infinite sets of inequalities which generalize all the known inequalities that can be used in the majorization step 
of the Approximating Hamiltonian method are derived. They provide upper bounds on the difference between 
the quadratic fluctuations of intensive observables of a JV-particle system and the corresponding Bogoliubov- 
Duhamel inner product. The novel feature is that, under sufficiently mild conditions, the upper bounds have 
the same form and order of magnitude with respect to N for all the quantities derived by a finite number of 
Mh' commutations of an original intensive observable with the Hamiltonian. The results are illustrated on two types 

of exactly solvable model systems: one with bounded separable attraction and the other containing interaction 
of a boson field with matter. 
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The Approximating Hamiltonian Method (AHM) provides a rigorous approach to the study of 
some classes of statistical mechanical systems in the thermodynamic limit. The method consists of 
the following interrelated ingredients: 

(i) Description of classes of model systems which admit a rigorous treatment in terms of a more 
simple approximating Hamiltonian; 

(ii) Rules according to which the approximating Hamiltonian is constructed from the original 
one; 

(iii) Mathematical techniques for derivation of bounds which prove the thermodynamic equiva- 
lence of the approximating and original Hamiltonians; 

(iv) Investigation of the thermodynamic and statistical properties of the system described by the 
approximating Hamiltonian. 

For the first time, the idea of the Approximating Hamiltonian method (AHM) has been sug- 
gested on a heuristic level by N.N. Bogoliubov in his paper on the theory of the weakly non-ideal 
Bose-gas He conjectured that, under the existence of Bose condensate in the system, the nor- 
malized (by the square root of the volume) creation/annihilation operators for bosons with zero 
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momentum can be replaced by c-numbers in the model Hamiltonian. The value of these complex- 
conjugate numbers was determined by thermodynamic arguments. Thus, the initial model Hamil- 
tonian was replaced by an approximating (trial) one, which has the advantage to be easily diago- 
nalized by a canonical u-v transformation. The further development of the AHM took place in the 
framework of the Bardeen-Cooper-Schrieffer (BCS) reduced model of superconductivity [2]. In 1957 
the corresponding approximating Hamiltonian was suggested and diagonalized by means of a u-v 
canonical transformation fs']. A little later, N.N. Bogoliubov [3| rigorously proved that the ground 
state energy density and the zero-temperature Green functions for the model and approximating 
Hamiltonians coincide in the thermodynamic limit. 

The foundations of the modern formulation of the AHM have been laid down by N.N. Bo- 
goliubov (Jr.), see and references therein. The essential generalization to classes of quantum 
systems with separable interaction, considered at nonzero temperatures, has been achieved by 
using the Bogoliubov variational principle for the free energy density and a special majorization 
technique based on integration over external sources. Some major restrictions on the applicability 
of the method, such as the quadratic form of the interaction Hamiltonian and its boundedness, 
were removed by further extensions of the AHM, see the review article Q and the books 0, @|. 

In this paper we generalize all the known inequalities which have been, or could have been, 
used in the majorization step (iii) of the AHM. They set different lower and upper bounds on 
the quadratic fluctuations, proportional to the difference between the original and approximating 
Hamiltonians, by terms involving new functionals of which the Bogoliubov-Duhamel inner product 
is a special case. Their place and role (not only in the AHM) is illustrated by the derivation of 
some new consequences pertaining to the infinitely coordinated anisotropic Heisenberg model and 
the Dicke superradiance model. These inequalities can be used for estimating the closeness of the 
Gibbs average values to the corresponding Bogoliubov-Duhamel inner product for a special class 
of observables. In some sense they are complementary to the Bogoliubov inequality [9] which has 
been used to exclude conventional superfluid and superconducting long-range order in one- and two- 
dimensional systems with gage invariant interactions. By exploiting it for setting upper bounds 
on the spontaneous magnetization or sublattice magnetization, Mermin and Wagner fio'l have 
rigorously proved the absence of long-range order in one- and two-dimensional isotropic Heisenberg 
models with finite-range interactions. Harris [11] has derived a lower bound on the symmetrized 
average value of the product of an operator A and its conjugate A\ which is a special case of 
the Bogoliubov inequality, but turns out to be sufficient for the derivation of the result of Mermin 
and Wagner. Alternative inequalities, setting a lower bound on the Bogoliubov-Duhamel inner 
product, have been used in the proof of the existence of spontaneous magnetization in a variety of 
quantum spin systems (l2| - [l5t . In the physical literature, the Bogoliubov-Duhamel inner product 
is also known as the Bogoliubov-Kubo-Mori scalar product. In addition, it plays an important role 
in the linear response theory [Tfij, Kondo problem [12], the so-called parameter estimation problem 
in quantum statistical mechanics and noncommutative probability theory see also |18l] and 
references therein. 



2. Systems with bounded separable attraction 

Here we consider a special class of quantum statistical models which are defined initially in 
a finite region A of the d-dimensional Euclidean space R'^ or the integer lattice Z"^. By ]A] we 
denote the volume of A in the first case, or the number of lattice sites in the latter case. The 
Hamiltonian Ha is defined as a self-adjoint operator in a separable Hilbert space H, and the 
corresponding free energy density /a[^a] is assumed to exist. For the sake of simplicity, we do 
not explicitly distinguish between a Hamiltonian Ha, describing a system with fixed number of 
particles TV in A, and the statistical operator Ha — M-A/" in the grand canonical ensemble, where ^ is 
the chemical potential and A/" is the number operator of particles. The density of the corresponding 
thermodynamic potential is given by 

/A[HA] = -(/31Al)-ilnZ[-«A], (2.1) 
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where Z[Hf^ is the partition function; in all cases the thermodynamic limit is denoted by t— lim. 
The norm of a bounded operator A is ||A||, the symbol remains for both the operator A and 
its adjoint A^ . As usual, [A, B] = AB — BA is the commutator of two operators. Average values in 
the Gibbs ensemble with the Hamiltonian T-L\ are defined as 

(■••)„, ^Tr(e-'^«-...)/^[HA]. (2.2) 

Let the Hamiltonian of the system in A be defined as a sum of two self-adjoint operators, 

nA=TA+UA, (2.3) 

where 7a = 7^ is a trace-class operator which generates the Gibbs semigroup {exp(— /37a)}^^o- In 
addition, we impose the condition that the density of the thermodynamic potential corresponding 
to 7a is bounded uniformly with respect to the volume of the system |A|, 

|/A[rA]KAfo. (2.4) 

A distinguishing feature of the models with bounded separable attraction is that the interaction 
Hamiltonian can be written as an extensive self-adjoint operator of the form Q 

n 

Ua = -\A\Y.9sAa,sAI,. (2.5) 

s=l 

Here > 0, s = 1, . . . ,n, are interaction parameters and the intensive observables Aa,s, s 
represent uniformly bounded local operators averaged over a region of the space (real space or 
conjugate momentum space). The uniform boundedness of these operators, 

II^a..||-PLIK^'^i> 5 = 1,...,", (2.6) 

where the constant Mi is independent of the volume |A|, is essential for the applicability of the 
method in the case under consideration. In the general framework one does not need explicit 
expressions for the operators 7a and Aa.s- It suffices to impose, in addition to (|2.4p and (|2.6p . the 
following general constraints (s, s' = 1, . . . , ri): 

||[^A,.,rA]|| Afi,T, (2.7) 

\\[Aa,s,A*J\\ ^ M2\A\-\ (2.8) 

The heuristic rule for construction of the approximating Hamiltonian H^^' (a) consists in lin- 
earization of the original interaction Hamiltonian with respect to the deviations of the intensive 
operators Aa^s, ^a s from some complex numbers: 

n 

T^f (a) - Ta - |A| 5s(as4,s + <^A.s - flsO, (2.9) 

5-1 

where a = (ai,...,a„) G C", and a* is the complex conjugate of Og. The latter numbers are 
considered to be variational parameters, chosen so as to minimize the contribution in the free 
energy density of the residual interaction Hamiltonian 

n 

7^W(a) ^Ha- (a) = -|A| ^^.(Aa,^ - a.)(Ai ^ - <) ^ 0. (2.10) 

The main result of the AHM for this class of model systems is summarized in the absolute 
minimum principle for the approximating free energy density: 

^ min/A[Hf (a)] - /a[Ha] < e(|A|), (2.11) 
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where e(|A|) as |A| oo. This estabhshes the thermodynamic equivalence of the free energy 
densities for the model Hamiltonian (|2.3p . (|2.5p . and the approximating one (12. 9p . 

The proof of (|2.1ip is carried out with the majorization technique developed by Bogoliubov Jr. 
It starts with the introduction of auxiliary external fields i/g and v* conjugate to the operators A]^ ^ 
and Aa,s, respectively, 

n n 

HaH ^7^A-|A|^(;.s<,+^^*^A,.), (a,^) (a)-|A|E('^^4,s+^r^A,s). (2.12) 

s=l s=l 

Next, lower and upper bounds on the difference in the free energy densities for the original and 
approximating Hamiltonians are set by the Bogoliubov variational principle and the inequalities 
following from it. In view of the non-positive definiteness of the residual interaction Hamiltonian 
Ha ''(a) = 'H\{v) — 'Hf\a, v), the application of the Bogoliubov inequalities yields the bounds 

^ min/A[Hf (a,z/)]-/A[HA(t^)] 

a 

71 

< E 9s{{Ak^s - (Aa,s)w^(.))(<s - (4,s>«a(.))>«a(.) . (2-13) 

s=l 

valid for all complex fields s — 1^2, . . .n. 

The fact that the absolute minimum of the free energy density /a ['Ha''' (a, v))] is attained at a 
finite value a = a\{v), which depends on the thermodynamic parameters of the Gibbs ensemble 
as well as on the size and shape of the domain A, has been proved by using the conditions (j2.4p 
and (EH). 

Here it is in place to mention that the attempt to prove directly that the correlation functions 
in the right-hand side of inequalities (|2.13p tends to zero as |A| — > oo may happen to be an 
impossible task. An efficient means for solving the problem provides the majorization technique 
of Bogoliubov (Jr.), see Instead of direct evaluation of the above correlation function, it uses 
inequalities which set upper bounds to averages of the form 

{SAA^sSA{ jnA{^) , SAa^s = Aa,s - (^a.s)«a(^) , (2.14) 

in terms involving second derivatives of the free energy density with respect to external fields lyg- 

^^MM^^^p\A\iSAA.s;SAA,s). (2.15) 

Here (A; B) denotes the Bogoliubov-Duhamel inner product on the algebra of observables A, B, . . . , 
defined as 

1 

{A; B)n = {Z[H]y^ J dr Tr [e-'^^^-^^^A^'e-'^^^B] . (2.16) 



In the remainder we consider the given observables A, B^C, . . . and fixed Hamiltonian H per- 
taining to a quantum system in a finite region A. Whenever no confusion could arise, for brevity 
of notation we will omit the subscripts H and A and the argument of the partition function Z . 

2.1. The Bogoluibov-Duhamel inner product 

The general properties of the functional (|2.16l) are considered in the book [31 and in the arti- 



cles |12l - ll4l . llTI . Il8l | . We warn the reader that some authors use definitions which differ from (I2.16P 
by a factor of /? and/or by involving the operator A instead of its adjoint A^. For our purposes 
it suffices to mention the following. The inner product [A; B) is conjugate symmetric, {A; B) = 
{B^;A^), antilinear (linear) in the first (second) argument, {A + aC;B) = {A;B) + a*{C;B) 
{{A; B + aC) — {A; B) -I- a{A; C)), it satisfies the relationship 

P{A-[H,B])u^{[A\B])u. (2.17) 
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The inner product (A; A) is nonnegative, {A; A) ^ 0, and convex, 

(l/2)(AAt +^t^). (2.18) 

Finally, in the case when either Aoi B (or both) commute with the Hamiltonian, the inner product 
{A; B) reduces to the bilinear complex-valued functional 

{A:B)n = {A^B)n. (2.19) 

Note that in some works the operator metric in the algebra of physical observables defined 
by (|2.19p is called Kubo-Martin-Schwinger metric, and the one defined by (I2.16P is called Bogoliubov- 
Kubo-Mori metric [ii,!!^. 

The further considerations are conveniently carried out by using the spectral representation of 
the Bogoliubov-Duhamel inner product. We assume that the Hamiltonian "H has a simple discrete 
spectrum only, {En, n — 1, 2, 3 ... } and denote the corresponding eigenfunctions by \n), i.e., T-Lln) — 
En\n}, n = 1,2,3.... By A^n = (to|A|7i) we denote the corresponding matrix element of an 
operator A. Then, the right-hand side of (j2.16p can be written as 

. P-/3-E™ _ P-/3B,. , 

{A;B)n = {ZA[H]r'J2'Al^B^n _ ^ . + (Za[H])-^ ^ c"^^"4;„i?„„ , (2.20) 

m,n ^ n 

where the prime in the double sum means that the term n = to is excluded. 
Our aim is to majorize the quadratic fluctuations 

{5AHA)^{A^A)u-\{A)\^ (2.21) 

by terms proportional to some power of the inner product 

{5 A- 5 A) = {A-A)~\{A)\^ 

„-PE^ _ -/3E„ 

= +(ZA[H])-^^e-^^"K„|^-|(A)r (2.22) 

m,n ^ ^ n 

Instead of (j2.2ip it is more convenient to consider the symmetrized form 

i(MtM + MM^) = ^{A^A + AA^) - \{A)\^. (2.23) 

Here we have 

^{A^A + AA^) = z-i^e-^^"i(|A™„|2-H|^„„|2) 

m,n 

= Z-i^'|^_pi(e-^^"+e-^^™) + ^-'El^"»l'e"^''"- (2.24) 

m,n n 

By comparing equations (|2.22p . (|2.23p and (I2.24p we obtain 

+ AA<) - iA;A) = Z-^Y.'\^,.,.f H- e-»~) - ^^j^^} ■ (2.25) 



Now, by using the identity 

^-PE^ + ^-0E„ ^ ^^-PE^ _ e-P^-)coih /^(-^"^ - -^") (2.26) 
we can express the difference (j2.25l) in two equivalent forms: 

-{A^A + AA'^) - {A- A) = V 'I A™ I' -577^ Er^(^™ cothX„„ - 1), (2.27) 



13003-5 



J.G. Brankov, N.S. Tonchev 



\{A^A + AA^) - {A;A) - ^ '|^_ri(e-'^^" + e"^^") fl - \ ) , (2.28) 



and 
1 

2^ / \ ' / I ""H 2^ ' \^ Xm,„cothX„ 

where 

X,nn^\li{E^-En). (2.29) 

The application of the elementary inequality a;cotha; ^ 1 to the right-hand side of (12.27^ 
immediately yields the convexity property (|2.18p . 

Different choices of the upper bound on the right-hand side of (|2.27p generate different inequal- 
ities. Thus, the inequality of Brooks Harris fll|, 

{A-A)^\{AA+ + A+A)^{A-A) + ^{[[A+,n],A]) (2.30) 



is obtained by setting 



1 < xcothx < 1 + ix^. (2.31) 
o 



On the other hand, if one uses another elementary inequality, 

1 < a;cotha; s: 1 + (2.32) 
and subsequently applies the Holder inequality, one obtains the result due to Ginibre (22|: 

{A- A) ^ ]^{AA++A+A) ^ {A- A) + \{{A- A)P{[[A+ A])}-^ . (2.33) 

A different choice of the parameters in the Holder inequality, followed by the implementation 
of the upper bound 

generates a symmetric version of the inequality due to Bogoliubov (Jr.) 

]^{AA+ + A+A) ^ [A- A) + i[(A; A)/3]2/3{([A+,H][H,A] + [H,A][A+,n])Y'\ (2.35) 

Due to relation (|2.15p . each of the above inequalities can be used in the AHM to majorize the 
quadratic fluctuations (j2.14l) in terms involving second derivatives of the free energy density with 
respect to the external fields Vs- However, note that in the zero temperature limit j3 oo the 
right-hand side of the simplest inequality (|2.30D diverges. 

3. Main inequalities 

To derive generalizations of the known inequalities involving the Bogoliubov-Duhamel inner 
product, a set of new functionals Fk{J; J), fc = 0, 1, 2, . . . , is defined by their spectral representation: 

Fk{J;J) :=Z-i^lJ™ip|e-^^' -(-l)V^^"|(/3l£;„-i?z|)'"^ (3.1) 

ml 

The specific choice of = 0, 1, 2, . . . is motivated by the relation of functionals (j3.ip to the Gibbs 
average values of some commutators and anticommutators involving the given operators J and 
Hamiltonian H. Indeed, 

(i) If fc = 2n, n = 0, 1, 2, 3, . . . , then 

ml 

= /32"(i?„;i?„) = /32"-i([i?+i?„_i-i?„_ii?+]), (3.2) 
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where, by definition, = Xj-u is a solution of the operator equation J — and 

Ro = Ro{J) = J, Ri=Ri{J)^[H,J], Rn = Rn{J)^[H,Rn-i{J)], n = 1,2,3,... 

(3.3) 

These observables have been introduced in [23j. 
(ii) If fc = 2n + 1, n = 0, 1, 2, 3, . . . , then 

F2n+i{J; J) = Mie-^""' + e-'5^™)[/3(i?„, - Ei)f^ 

7nl 

= /32"([i?„i?+ + i?+i?„]). (3.4) 

In particular, 

Fo(J;J) = (J;J), Fi{J;J) = {JJ+ + J+J), F2(J; J) = f3{[[J+ ,H], J]), 
F3{J;J)^(5^[J+,H][n,J] + [H,J][J+,H]). (3.5) 

The functionals p.ip will be used to generalize all the known inequalities used in the AHM. 

3.1. Generalization of the l-iarris inequality 

By using the identity (|2.26p . one can rewrite the equality (|3.4p in the form 

/3(S„-S0, 

\'-iml\ ■ — e ; cuui 

ml 

Now, from the elementary inequalities (j2.3ip it follows that 



f2n+i(J; J) ^ iJrnil'ie-^''' - e"^^-) coth [/?(£;,„ - g,)]^"- (3-6) 



F2n{J; J) ^ ^F2n+l{J; J) ^ F2n{J \ J) + ^i^2n+2(J; J)- (3.7) 

This is a generalization of the Brooks Harris inequality (j2.30p , since the latter is recovered in 
the particular case of n = 0. 

3.2. Generalization of the Plechko inequalities 

The application of the elementary inequalities (|2.32p to the right-hand side of p.6p leads to 

F2n{J;J) < i^^2„+i(J;J) s=;F2„(J;J) + (2Z)-i^|J„,zne-^^' -e-^^-|[/3(i?,„-i?0]'". (3.8) 

ral 

The problem now is that, due to the absolute value of the difference of the two Gibbs exponents, 
there is no apparent interpretation of the sum in the right-hand side of (|3.8p in terms of average 
values. A known way to overcome this difficulty is based on the application of the Holder inequality 

C\ i/p / \ 1/9 

By setting k = (m, /) and 



|2 



e ' — e 



/3{E^ - El) 



^ ^ lj{E„,-Ei) \P\Era-Ei)\r-+'^'' \ , (3.10) 
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we obtain 

ml 

One of the possible choices of p and q here, namely even integer q — 2k (hence, p ~ 2k/(2k— 1)) 
leads to the set of generalized Ginibre inequalities (fc = 1, 2, 3, . . . ): 

F2niJ;J) ^ iF2„+i(J;J) ^F2„(J;J) + i(J;J)(2'=-i)/2'=[F2fc(2„+i)(J; J)]!/^^ (3.12) 
Atn = the above set reduces to a symmetric version of the inequalities obtained by Plechko (2l| : 
(J; J) < i(JJ+ + J+ J) ^ (J; J) + i(J; J)(2'=-i)/2fe/3(i?,; i?,)i/2fe, (fc ^ l, 2, 3, . . . )■ (3.13) 
Hence, in the particular case of fc = 1 one obtains the Ginibre inequality (j2.33p . 

3.3. Generalization of the Bogoliubov (Jr.)-Plechlco-Repnikov inequalities 

If in (I3.1ip one chooses odd g = 2fc + 1, hence, p — {2k + I) /2k, then for the sum in the 
right-hand side one can use the upper bound 

ml '■^ 

^ |J™n'|e"^^' +e-^^"|[/3|i?™ = F2i2nk+n+k)+i{J;J)- (3.14) 

7nl 

Thus one obtains the set of mequahties (/c = 1, 2, 3, . . . ) 

ii^2„+i(J; J) F2n{J;J) + i(J; J)''=/(''=+^)[F2(2„fc+„+fe)+i(J; J)]i/(2'=+i). (3.15) 

At n = these reduce to a symmetric version of the set of inequalities obtained by 
Bogoliubov Jr., Plechko and Repnikov (23| : 

^{JJ+ + J+J) <: (J;J) + i(J; J)2fc/(2'=+i){/32fc(i?,i?+ + i?+i^,)}i/(2fc+i). (3.16) 

The symmetric version of the Bogoliubov's (Jr.) inequality (j2.35p follows from here in the 
particular case of A: = 1. 

3.4. Alternative sets of inequalities 

By applying the Holder inequality p.9p to the second term in the right-hand side of p.7p under 
the substitution 

\ymi\ = {iJmil'le-^''^ -e-f^^-imrn-mir+^/^y^', (3.17) 
instead of p.lOp . one can in parallel derive two new sets of inequalities. Thus we obtain first 

(2Z)-i ^ I J^zHe-^^' - e-^^™|[/?(£;„, - 

ml 

^ ^[^2„(J; J)]^/^ \Z-'Y.\J^^\" %(^E^ 1 E,) mEm-EAr+''(^-'A . (3.18) 
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Now there are two choices of p, one of which yields l/{p — 1) odd integer, and the other - even 
integer. In the first case we set — 1) = 2fc — 1, fc = 1, 2, 3, . . . , which imphes p = 2k /{2k — 1), 
q — 2k. Then, from p.lSp and p.7p the foUowing set of inequahties follows 

iF2„+i(J; J) «CF2„(J; J) + i[F2„(J;J)](2'=-l)/2'=[F2(„+fe)(J;J)]l/2^ (fc = 1, 2, 3, . . . )■ (3.19) 

In the second case we set l/{p — 1) = 2k, k = 1,2,3,..., which implies p = {2k + l)/2k, 
g = 2fc + 1 . Then we can use the upper bound 

'ml 

^ Z-i^ |J™n'|e-^^' +e-'3^-l[/3|i?™ - = i^2(„+fc)+i(J; J), (3.20) 

ml 

to obtain another set of inequalities 

iF2n+i(J; J) < F2n{J; J) + i[F2„( J; J)]2'=/(2fe+i) [F2(„+,)+i ( J; J)]i/(2'=+i), {k = 1, 2, 3, . . . ). 

(3.21) 

Note that (|3.19p is different from (|3.12p but at n = it reduces to the Plechko inequali- 
ties (I3.13p . Therefore, this set of inequalities can also be considered as a generalization of the 
Ginibre inequality (|2.33p . 

Similarly, for general n, equation p.2ip differs from (|3.15p but reduces to the Bogoliubov (Jr.)- 
Plechko-Repnikov inequalities (I3.16P at n = 0. Hence, the set of inequalities p.2ip is also a gener- 
alization of the Bogoliubov (Jr.) inequality (|2.35p . 

A comment is in order here. Due to the property Rn+k{J) = Rk{Rn{J)), in terms of the 
operators Bn = Rn{J) one has 

F2„(J; J) = /32"(B„;B„), ^^2„+i(J; J) = /?'"(S„B+ + B+B„), 

F2(n+k) + l{J\J) - (i?fe(B„)4(^n) +4(^n)^'o(^n))- (3-22) 

Therefore, inequalities p.7p take exactly the form of the inequality of Brooks Harris (|2.30p with 
A replaced by i3„, inequalities (|3.19p . respectively (|3.2ip . take exactly the form of a symmet- 
ric version of the Plechko inequalities (j3.13l) . respectively, the Bogoliubov (Jr.)-Plechko-Repnikov 
inequalities (I3.16p . with J replaced by i3„. 

Notably, under the above substitution, the generalized Ginibre inequalities (|3.12p and the gen- 
eralized Bogoliubov (Jr.)-Plechko-Repnikov inequalities p.lSp do not reduce to any of the known 
types of inequalities, except in the particular case of n = 0. 

3.5. General features and comparison of upper bounds 

Obviously, the main inequalities have been derived under different approximations. The most 
direct is the derivation of generalized Harris inequalities - it is based upon the single elementary 
upper bound (|2.31l) . Next, the generalized Ginibre inequalities are derived by first using the el- 
ementary upper bound (|2.32p . followed by the application of the Holder inequality (|3.9p under 
a special choice of the parameters: p — 2k/ {2k — 1) and q — 2k. An alternative choice of these 
parameters, p = {2k + l)/2k and q — 2k+\, requires the use of an additional, rather crude upper 
bound (|2.34p . in order to derive the Bogoliubov (Jr.)-Plechko-Repnikov inequalities (|3.15p . 

The characteristic feature of our generalized inequalities is that, under sufficiently mild condi- 
tions on the Hamiltonian of the A^-particle system and the bounded intensive observable J^v, 
they provide the upper bounds on the non-negative difference 

A„(JAr) = ]^{R+{JN)Rn{JN) + Rn{JN)RWN))nM - (i?n( ^iv); i?n( J7v))«„ (3.23) 
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of equal form and equal order of magnitude (with respect to N) for all finite n — 0,1,2,..., 
i.e., for any finite set of observables Jn, Ri{J]s[) = [Hn, Jn], R2{Jn) = \HntRi{Jn)\, ■■ ■ ■ Note 
that the left-hand side inequality in (j3.23l) is a generalization of the convexity property of the 
Bogoliubov-Duhamel inner product (|2.18p . 

The required conditions are (n = 1, 2, 3, ... ): 

\{Jn)hA < 

F2n{JN]JN) = /3'"-^([i?+(JAr),i?„_l(JAr)])«„ S;0(A^-1), 
F2n+l{JN]JN) = l3^''{Rt{JN)Rn{JN)+Rn{JN)Rt{JN))'H^ ^0{l). (3.24) 

They are generally satisfied for extensive Hamiltonians H n with bounded interaction and bounded 
intensive observables Jn, which are arithmetic averages over the particles of the system of some 
local observables. An example will be given in the end of this section. 

The above conditions may also hold in some cases when Ti^ and i?„(JAr), n = 0,1,2,..., 
contain unbounded operators. As an example in section 14.21 we will consider the Dicke model of 
superradiance [25!|, for which the proper unbounded counterpart of Jjy will be found. 

Under conditions p.24p . our generalized inequalities yield the following upper bounds for all 
finite n = 0, 1, 2, . . . : 

1. Generalized Harris inequality (|3.7p 

A„(Jw) ^ 0{N-'). (3.25) 

2. Generalized Plechko inequalities p.l2p 

A„(Jiv) ^ (JAr;Jjv)^"=-''/'"0(^-'/''=), (fc- 1,2,3,...), (3.26) 
which at fc = 1 reduce to the Ginibre inequality 

A„( Jat) ^ {Jn; JnY'^0{N-^'^). (3.27) 

3. Generalized Bogoliubov (Jr.)-Plechko-Repnikov inequalities (|3.15p 

A„(Ja') (Jiv;Jjv)"=/^''+'^0(l), (fc = 1,2,3,...), (3.28) 
which at /c = 1 reduce to the Bogoliubov (Jr.) inequality 

A„(Jw) < (Jat; Jw)2/3o(i). (3.29) 

Due to the relationship between the Bogoliubov-Duhamel inner product and the susceptibility 
of the system with respect to the external field conjugate to the observable Jn, see (j2.15p . we 
can compare the different upper bounds in the region of parameters in which the susceptibility is 
bounded. We see that the generalized Harris, Plechko and Ginibre inequalities yield upper bounds 
of the order 0(iV~^), while the generalized Bogoliubov (Jr.)-Plechko-Repnikov inequalities yield 
upper bounds of the order 0(Ar-2fc/(2fc+i))^ fc = 1, 2, 3, . . . . 

Finally, to illustrate the validity of conditions (|3.24p and the explicit form of the first generalized 
observables i?„(Jjv), we consider a simple model with separable attraction built upon bounded 
operators. Let the model system contain N spins and have the Hamiltonian (normalized by fceT) 

PUn - -Ng^ [Jfj f - Ngy {J^f - A^h • Jat , (3.30) 

where g^, gy > are dimensionless coupling constants, h = {hx, hy, hz} is the vector of the external 
magnetic field, J^r = {J^, J^, J^} is the vector operator of the average spin with uniformly 
bounded in N components J^, 

1 ^ 

•^^^TV^"^*"' Oi^x,y,z, (3.31) 

i=l 
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where cr" are the standard PauU matrices. 

By using the commutation relations for the spin operators, we obtain the observables 

R,{Jff) = [mN, J^] = ngviJlJN + JnJI) + hyjf, - KJl], (3.32) 

and 

i?2(J^) = [mN. Ri{Jf,)] = 45,(5, - gMJlJN + '^JIJn'JI + J%JlJl) 

+ 4(5y - gx)hy{J%Jl + J^) - 45,/i,(Jj^ + J^) 
which, by definition, have zero average values. Next 



+ A{hl + hl)J^ ~ AKhyjy, - 4h^h,r^ , (3.33) 



F,iJ%;r^) ^ (i?i(J^);i?i(Jj^))==([i?+(J^),J^]) 

= ^{29y[{iJ'Nr)-{iJNn+hy{Jf,) + h.{J^)}, (3.34) 



and 



FM;r^) ^ {[Rt{r^)RM) + R,iJmtiJN)]) 

= {[gyiJ^J' + JNJN) + hyJN"h,JN?). (3.35) 
Rather lengthy but straightforward calculations show that 

|F4(J^; J^)| = \{[RtiJN),Ri{J%m ^ 0{N-^). (3.36) 

One can readily extend the above results and show that for all finite n the following inequalities 
hold: Ki?„(J^))| ^ \\RniJ^)\\ ^ 0(1), |F2„(J^; J^)| ^ 0(7V-i), and |F2„+i(J^; J^)| 0(1). 
Thus, conditions p.24p and, hence, the upper bounds (l)-(3) are fulfilled for the considered model. 

4. Systems of matter interacting with Boson fields 

4.1. Models and their treatment by the AHM 

Here are two examples of models in solid state physics, which belong to this class: (a) The 
Dicke model of superradiance [25], solved exactly in [26, 27 1, and by the AHM in ^2^. The model 



has been generalized to include interactions with both electromagnetic field and phonons [.29, 30| , 
and to the case of infinitely many modes of the electromagnetic field [U . A recent review of the 
thermodynamic properties of the original Dicke model and its generalizations is given in [3^ . (b) 
The Mattis-Langer model of structural instability [33] , solved exactly by the AHM in [sH . A class 
of models, including as a special case the Dicke model, has been considered in the framework of 
the AHM by Bogoliubov (Jr.) and Plechko [35]. The one-dimensional case with countably infinite 
set of phonon modes has been solved by means of theta- function integration in [36]. 

The model Hamiltonian is defined on the tensor product of two Hilbert spaces, one for the 
subsystem describing matter (e.g., electrons in a solid, spins on a lattice), and the other for the 
Boson field (lattice vibrations, electromagnetic field). In the second quantization representation, 
the creation, 6|, and annihilation, feg, operators of the Boson field modes (labeled by the subscript s) 
satisfy the canonical commutation relations 

[h,bj]^6,.i, [b,,bi] = [blbl]=0, (4.f) 

for all the allowed s and I. Since the Boson operators are unbounded, it is not possible to obtain 
easy bounds on their average values in terms of Hilbert-space norm. 
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Another characteristic feature of this class of models is that exact solvability by the AHM is 
possible when the interaction with only a finite (or growing slower than the volume V , as V ^ oo) 
number of Boson modes is taken into account. The typical model Hamiltonian has the form 

n n 

Ha = Ta + ^ b% + As(&sA^ + blAs,A)- (4.2) 

s=l s=l 

Here the operators 7a = and As.a, s — 1, . . . ,ti, refer to the matter subsystem and satisfy the 
general conditions (|2.4p and (|2.6p . The second term in the right-hand side of (14. 2p describes a finite 
number of free Boson modes, s = 1, . . . , n, in the space domain A of volume V. For the sake of 
simplicity, the energies Wg > and the interaction constants Ag G i? are taken to be independent 
of the volume |A|. The mathematical definition of the above Hamiltonian is given in |37|. 

The corresponding approximating Hamiltonian depends on a set of complex numbers r] = 
(?7i, . . . , rjn) G C" and has the form, see .28] . 



Hf\v) - ^c., b% + Ta -Vj2i^l/u:,){v,Al j, + <^,a - <^s), (4.3) 



,io) 

s=l s=l 

where 



11 = bl + V^/^\t , 6. = 6. + V^'^\, , (4.4) 

are the creation/annihilation operators for a subsystem of free shifted bosons. The application 
of the Bogoliubov inequalities to the difference of the free energy densities for the model and 
approximating Hamiltonians yields 

mTm-mK]^Q. (4.5) 

for all ?/ G C". Therefore, the best approximation is reached at r/ = f/A, where t/a satisfies the 
absolute minimum condition 

/[^f (^A)]=min/[Hf (ry)]. (4.6) 
Note that the free energy density of the free bosons. 



/ 



b% 



^i:in(l-e-^-), (4.7) 



is independent of the parameters 77 = (771, . . . , ?]„) G C" and vanishes in the thermodynamic limit 
as 0{V-^). 

The proof of the thermodynamic equivalence of the free energy densities J[H!"^\f\A)\ and /["Ha] 
goes again through the introduction of sources of the boson fields, 

n 

HA(t^) = HA-V^'^Y.^v*h, + v,bl), 

n 

nf\il,i^) = (r;)-W/2^(^r6. + z.g&t), (4.8) 

s=l 

where v = (i/i, . . . , z/„) G C". 

Now, from the Bogoliubov inequalities and a subsequent use of an elementary upper bound, 
one obtains 

n 

«C m\nf[nf{ii,v)] - mA{i^)] ^ A3(56sM,,A + ^61M,,a)«^H 

S — 1 

n , 2 ^ 

LJc 

s=l ^ s=l 
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where 

6As^A = As^A- {As,a)-Ha{i^), Sbs^bs- {bs)-HA{,,), (4-10) 
and 7 £ (0, 1) is a free parameter. 

Due to the boundedness of the operators As^a, AI ^, see conditions (|2.6p . the first sum in the 
right-hand side of the last inequality (|4.9p is bounded from above by 

V^-i/2-7„Af2^ax(A2/ws). (4.11) 

s 

The quadratic fluctuations of the boson fields, (SblSbs) -UAi^)-! ^re to be majorized by terms 
proportional to powers of the Bogoliubov-Duhamel inner product 

By using the Ginibre inequality (|2.33p and choosing 7 = 1/3 in (|4.9p . the subsequent applica- 
tion of the majorization technique due to Bogoliubov Jr. with the following bounds on the first 
derivatives (s = 1, . . . , n) 



= 1^-1/2 1 (6f = u:-'\\Mtiv)nAiu) - vf\ < w-i(AsMi + \vf\), (4.13) 



dv* 

made it possible to prove that 

\m\nmf\Ti)]-m^]\^^v, (4.14) 

where ey = 0(y~^/'^) ^ as y ^ 00. Under additional conditions on the double commutators 
between the different {^f^ : s — 1, . . . ,n} and on the commutator of i?i(A^^) with ^, the 
above estimate was improved up to ev = Oiy-^''') 0. 

4.2. Application of the generalized inequalities to the Dicke model 

In the remainder, by using the equality (l3. [23|: 

P{X;[n,B])n^{[X\B])n, (4.15) 

we shall derive a variety of explicit relationships between Bogoliubov-Duhamel inner products 
and usual thermal averages for different observables of the basic single-mode Dicke model in the 
rotating wave approximation. In the latter case, the Hamiltonian has the form (|4.2p with n = 1, 



N N N 

i=i i=i i=i 



where crj^ = ■^{<Jj ± icrj) and cr| are the Pauli matrices. From the above definitions the following 
commutation relations follow 

[r,A]=eA, [A^,A]^^T. (4.17) 

By direct computation we obtain 

Ri{b) EE [H, b] = -(a;6 + V^^'^XA), (4.18) 

and, since {[H,b])-H — 0, we obtain a well known equality between average values of observables 
pertaining to different subsystems (boson field and matter): 

{b)^ = -V'/^-{A)n. (4.19) 

OJ 
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Therefore, in the case of non- vanishing polarization in the matter subsystem, (A)^ ^ 0, the average 
value of the boson annihilation (as well as creation) operator will behave as = 0(y^/^). Hence, 
we deduce that the unbounded counterpart of the observable Jat in conditions (|3.24p should be the 
normalized operator V^^l'^h. Keeping this in mind, we further calculate 



2A2 



(4.20) 



and 



F^iy-^'H-v-^'H) 



(4.21) 



V-\b% + bb^)-H + V-^^^-{b^A + bA^)n + K{A^A + AA^)h 



Fi{V-^/\V-^'%) = P^[Rl{V-^/%),Ri{V-^'%)])-H 



V 



2A^ / e \ , , 2A3 , , 



(4.22) 



(4.23) 



The right-hand side of these expressions can be evaluated using some relationships between 
average values of different observables which follow from (j4.15p . Thus, by setting B = b and 
X = b, we obtain, 

{b;b)n^^-V'^''-{b;A)n. (4.24) 

PUJ U! 

Next, from B = b and X = A it follows that 



-{A;b)n^V'^'-iA;A)n. 
Under the alternative choice B = and X = m. (|4.15p . we derive 

{b^-M)u^^-V'''-{b^;A^)n. 
which, due to the conjugate symmetry {A; B) — {B'^;A^)^ is equivalent to 



{b;b)n^^-V'/'^{A; 



)H ■ 



(4.25) 



(4.26) 



(4.27) 



By comparing this equality with (|4.24p we conclude that (6; A)-^ = {A;b)-u- Then, taking into 
account (j4.25p we derive the important relation 



{b;b)n^^ + V^{A;A)n.. 

puj CJ^ 



(4.28) 



which was obtained in [32] using gage invariance arguments. 

Proceeding further with the evaluation of (|4.22l) and (I4.23p . we note that since the right-hand 
side of (j4.23p is real, one must have 

{b'^A)n = {bA^)H ■ (4.29) 
Therefore, the application of the Schwarz inequality and the relation (|4.28p yield the estimate 

\{V-^'%'^ A)u\ = {V-^/HA^)h\ < {V-^b%)]{^\{AA^)]{''\ < 0(1). (4.30) 
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Finally, the application of the Harris inequality (|2.30|) in terms of b gives 

(6; b)n ^ {h%)H + \^ip- h)n + ^/3a; , (4.31) 

which, in view of (|4.28p . implies {h^h)u = 0(V"), provided that {A\ A)-u = 0(1). Taking into account 
Pll = 0(1) and ||r|| = 0(V), we conclude that F2iV-^/^b;V~^/^b) and F4{V-^/%;V-^/^b) are 
bounded from above by 0{V-^), while 1^^3(^-1/25.^-1/25)1 <^ Qf^^y 

Thus, we have proved that V'^^^b and FniV^'^^'^b^V'^/'^b) for n = 2,3,4 satisfy condi- 
tions ([X^ . Therefore, bounds of the form (1X2^ . (IX?7)) . and (IX^ hold true at n = 1 even 
for the unbounded operator Jy = V~^^^b. 



5. Discussion 

The AHM makes it possible to rigorously obtain the exact thermodynamic properties of diverse 
classes of model systems in quantum statistical mechanics. Each class of models requires a certain 
structure of the interaction Hamiltonian and special properties of the operators in terms of which its 
structure is defined. Initially, the mathematical technique for derivation of bounds which prove the 
thermodynamic equivalence of the approximating and original Hamiltonians has been developed for 
systems with interaction Hamiltonians constructed with bounded operators (see section 2). Later, 
it was extended to the case of unbounded Bose operators (see section 4) . This technique essentially 
exploits the possibility of estimating the correlation functions in (|2.14p or (|4.9p from above by ex- 
pressions containing the Bogoliubov-Duhamel inner product in combination with average values of 
certain operators. As a rule, the latter are estimated by norm, which requires conditions (|2.4l) - (|2.8p . 
It is well known that this procedure is not unique 0, @l • The distinction comes from the different 
inequalities used. At first, it was Bogoliubov's Jr. inequality that was used in the case of bounded 
operators. Latter it was realized that Ginibre's inequality (|2.33p is more convenient to be used, es- 
pecially in the case of interactions involving Bose operators [7|, see also section 4. Moreover, it made 
it possible to derive a better estimate, namely 0{\A\~^^^), in both cases of bounded and unbounded 
operators However, in the latter case, two additional sufficient conditions have been imposed. 

The interest that stimulates the search and use of different inequalities was prompted by the 
wish to improve the upper bounds on the difference of the model and approximating free en ergy 
densities and/or to enlarge the class of model systems rigorously solved by the AHM 12 ll |23|. 



In section 3, we have derived generalizations of the inequalities introduced in [21| and [23|. The 
novel point is that a set of new functionals Fk{J; J), = 0, 1, 2, . . . , were defined by their spectral 
representation p.ip . This made it possible to obtain different upper bounds in a unified fashion. 

Beyond the use in the majorization step of the AHM, our new upper bounds p.25p - (l3.29p on 
the difference between the Gibbs average values of a class of observables of a many-particle system 
and the corresponding Bogoliubov-Duhamel inner products may find wider applications. Under 
sufficiently mild conditions, these upper bounds have the same form and order of magnitude with 
respect to the number of particles (or volume) for all the quantities derived by a finite number of 
commutations i?„( J) of an original intensive observable J with the Hamiltonian of the system Ti. 

In addition, we have obtained important relationship between average values of the different 
observables in the framework of two types of exactly solved by the AHM model systems: one with 
bounded separable attraction - the infinitely coordinated anisotropic Heisenberg model, and the 
other, the Dicke superradiance model, describing interaction of a boson field with a subsystem of 
matter. 
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Vsara/ibHeHi HepiBHOCji fl/ia BHyrpiiuHboro flo5yTKy 
BorojiK)6oBa-flK)aMe/i5i is sacTOcysaHHJiMM b Merofli 
anpoKCMMyKDHoro raivii/ibTOHiaHy 

fl>K. r. BpaHKO^, H.C. ToHHeP 

/la6opaTopiiR xeopeTWHHoT 4>l3i/iKi/i iM. BorojiKiBoBa, 05'eflHaHi/iM iHCTHTyx aflepHnx floc/iifl>KeHb, 
141980 fly6Ha, PoclKcbKa tpeflepai4ia 

iHCTMxyT MexaHiKW, BoyirapcbKa aKafleMin HayK, Cocjjia, Bo/irapia 

iHCTi/ixyT 4>i3WKM TBepflOfo xi/ia, BojirapcbKa aKafleMia HayK, Co4)in, Bo/irapi^i 

OrpuiMaHO HecKiHHeHHi Ha6opH HeplBHOcreti, sikI y3ara;ibHHDHDTb Bci BiflOMi HepiBHOcri, mo MO>KyTb 6yTH 
BMKopwcxaHi Ha erani Ma>KopyBaHH5i Mejofly anpoKci/iMyKDHoro rawijibTOHiaHy Bohw 3a6e3neHyK)Tb sepxHi 

rpaHWL4i Ha pi3HMLlK) Mi>K KBaflpaTUHHI/lMH cjjiiy KTya L4i51 Ml/1 ihTeHCI/IBHUX CnOCTepe>KyBaHWX W-HaCTWHKOBoT cw- 

CTeMM i BiflnoBiflHoro BHyTpiiiiHboro flo6yTKy Boro/iio6oBa-flioaMe/i5i. Hoboio pucoK) e xe, lu,o npn flocxaTHbo 
M'iRKHx yMOBax Bepxni rpaHHL4i MaioTb oflhaKOBy 4>opMy i nopjiflOK BejinHUHH no BiflHoiueHHKi flo N fl/ia Bcix 
BejiMHMH, oxpuMaHMX ujjiaxoM CKiHHeHHoro HHC/ia nepecTanoBOK noHajKOBoT iHTencuBoT cnocTepe>KyBaHoY 3 
raMijibTOHianoM. Pe3y/ibTaTW i/iKDCTpyioTbC5i na flBox Ti/inax tohho po3B'a3yBaHMX MOflejiew: oflnieT 3 o6Me- 
meH\AM cenapa6e/ibHHM npuxaraHMJiM ja iHLUoT, L140 MicTi/iTb B3aeMOfliK) 5o30HHoro no/ia 3 MaxepieKi. 

K/iiOHOBi c/ioea: KopejifmiUHi (pyHKLiiY, BHyTpiwH'm floByroK EoroMio6oBa-flioaMe/ifi, HepiBHOcri 

CT3TI/ICTI/l'-IHOy MeXaHIKtl, MeTOfl anpOKCHMyHDHOFO raMijIbTOHiaHy, TOHHO po3B'ft3yBaHi MOfteni 
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